Ferromagnetism in one-dimensional metals: Breakdown of the Hartree-Fock 
approximation and possible first-order phase transition 
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We calculate the Gibbs potential r(M) of a one-dimensional metal at constant magnetization 
M to second order in the screened electron-electron interaction U. At zero temperature we find 
that r(A'f) contains non-analytic corrections proportional to Af^ln|Af| and \M\^, implying that a 
possible paramagnetic-ferromagnetic quantum phase transition in one-dimensional metals must be 
first order. 
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I. INTRODUCTION 



Can the ground state of a one-dimensional (ID) clean 
metal exhibit spontaneous ferromagnetism? A rigorous 
theorem due to Lieb and Mattis^ impUes that the an- 
swer to this question is "no" for ID continuum models as 
well as for ID lattice models with nearest neighbor hop- 
ping and interactions involving only densities. However, 
the Lieb-Mattis theorem docs not apply to lattice models 
with longer range hoppings. Indeed, some time ago Daul 
and Noack^ presented numerical evidence that the ID 
Hubbard model with nearest- and next-nearest neighbor 
hopping has a ferromagnetic ground state in a substantial 
range of densities and on-site interactions U . Given the 
stability of the ferromagnetic ground state in II? in a cer- 
tain parameter regime, one might want to know the crit- 
ical behavior of the system close to the quantum phase 
transition separating the paramagnetic from the ferro- 
magnetic regime. Some physical properties of ferromag- 
netic metals in 11? have been studied in several recent 
works^ii^. However, the fluctuation corrections to the 
Hartree-Fock approximation for the ground state energy 
have not been thoroughly investigated. In view of the fact 
that the metallic state in ID is a Luttinger liquid, it is not 
clear whether the Hartree-Fock scenario of a second order 
phase transition to a ferromagnetic state for sufficiently 
strong interaction is at least qualitatively correct. In this 
work we shall therefore calculate the leading (second or- 
der in U) correction to the Hartree-Fock appoximation 
for the Gibbs-potential r(M) at constant magnetization 
M. We find non- analytic terms which completely in- 
validate the Hartree-Fock prediction and imply that, if 
the paramagnetic-ferromagnetic quantum phase transi- 
tion exists in a ID metal, then it must be first order. 
Recently Belitz and Kirkpatrick^ came to a similar con- 
clusion about the paramagnetic-ferromagnetic transition 
in clean itinerant ferromagnets in dimensions 1 < Z? < 3. 

Although most of our considerations are fairly general 
and independent of any specific model, for explicit cal- 
culations we shall use the Hubbard model with nearest 
neighbor hopping t and next-nearest neighbor hopping t' 
on a one-dimensional lattice with = L/a sites and 



lattice spacing a. The hamiltonian is 



(1.1) 



k{7 



Here Cka = Nl ''''^'Cia- annihilates an electron 

with momentum k and spin cr, the operator fiio- = cJ^Cio- 
counts the number of electrons with spin a at lattice site 
Xi, and the energy dispersion is 



efe = — 2icos(fca) — 2i'cos(2fca) 



(1.2) 



II. RENORMALIZED PERTURBATION 
THEORY AT CONSTANT MAGNETIZATION 

A. General considerations 

In the presence of a uniform magnetic field h the grand 
canonical potential is 

n{ii,h) = -rin|Trexp[-/3(i? - fiN - hM)]^ , (2.1) 

where /i is the chemical potential, T — 1//3 is the tem- 
perature, N = X^fco- ^l:a'^k<j is the particle number oper- 
ator, and the operator M = J2ka ^'^ka'^^'^ represents the 
uniform magnetization. The expectation values of these 
operators are 

N = {N) = -dn/dfi , M = (M) = -dn/dh . (2.2) 

To study spontaneous symmetry breaking, it is more 
convenient to work with the corresponding Gibbs 
potentialiSiS 

r(iV, M) = n{fi{N, M), h{N, M)) +fj,N + hM , (2.3) 

which is a function of N and M. Here iJ,{N,M) and 
h{N,M) should be calculated by inverting Eqs. (|2.2|l . 
Perturbation theory generates an expansion of r(7V, M) 
in powers of the relevant dimensionless interaction, the 
so-called Stoner factor^" 



(2.4) 



2 



where Vq = = 0) is the density of states (DOS) per 
spin projection at the Fermi energy of the non-interacting 
system in the absence of a magnetic field. We define the 
energy-dependent DOS via 

Ho^^y.m-^k+fi), (2.5) 



N, 



k 



which has units of inverse energy. It is convenient to 
measure particle number and magnetization in reduced 
units 

(2.6) 



2Nr 



2Nr 



so that < n < 1 and — ^ < m < i. Note that the 
physical range of m- values is \m\ < n for n < 1/2 and 
\m\ < 1 — n for n > 1/2, because the magnetization M 
cannot exceed the number TV of electrons for a less than 
half-filled band or the number 2Nl — N oi holes for a 
more than half filled band. Depending on the values of 
the two dimensionless parameters n and 



7 = 4t7i , 



(2.7) 



the Fermi surface of our model without interaction and 
magnetic field consists of two or four discrete points. The 
different regimes are shown in Fig. ^ (A similar figure 
can be found in Ref. 0.) For I7I < 1 the Fermi surface 
has two Fermi points for all fillings < n < 1. For 
7 > 1, i.e. for t' > t/A, the Fermi surface has two points 
for < rt < nc(7), and four points for nc(7) < 71 < 1. 
The critical filling separating these regimes is 



/ N 1 (2 

rLcy^) ~ — arccos 1 

TT V 7 



7 > 1 



(2.8) 



On the other hand, for 7 < — 1, corresponding to t' < 
—t/A, the Fermi surface has four points at low fillings 
< 71 < He (7), and two points at larger fillings nc(7) < 
n < 1, where now 



ndl) = - arccos (7^ 
TT VI7I 



1 



7 < -1 



(2.9) 



Besides the usual singularities for n — > and n ^ 1, for 
I7I > 1 the DOS at the Fermi energy exhibits additional 
one-sided singularities at the critical fillings nc(7), which 
are related to vanishing Fermi velocities when the Fermi 
surface topology changes discontinuously. 

In the regime where the Fermi surface consists of two 
points the DOS at the Fermi energy is 



1 



TTVp 27ri sin(7rn)[l -|- 7COs(7rn)] 



(2.10) 



where is the Fermi velocity. The DOS is more com- 
plicated in the regime where the Fermi surface has four 
points. For 7 > 1 and tt. > ric we find 



2 



1 



tan(7r(l - n)/2) 
1 -I- cos(7r(l — 



n)) 



[l + cos(7r(l-n))]2-4/72 



(2.11) 




FIG. 1: Topology of the Fermi surface for the energy disper- 
sion l|1.2^ as a function of 7 = 4t'/i and filling n. The Fermi 
surface consists of two points in the shaded regime and has 
four points in the white regime. The insets show the typical 
dispersions with the position of the chemical potential as a 
horizontal line. 



and for 7 < — 1 and n < ric 



Vq = 



1 



tan(7rn/2) 
1 -|- cos(7m) 



[1 + cos(7rn)]2 - 4/72 



(2.12) 



Graphs of the DOS at the Fermi energy for different val- 
ues of 7 are shown in Fig. |21 For simplicity we assume 
in this work that in the absence of a magnetic field and 
ferromagnetic symmetry breaking the true Fermi surface 
consists only two points kp and —kp^ each of which is 
split by the magnetic field h (or for finite to) into two 
spin-dependent points ka and —ka^ where a =t7j,- The 
calculations in this work are therefore restricted to the 
shaded regime shown in Fig.^ We emphasize that ka are 
the true Fermi momenta of the interacting system, which 
have to be determined self-consistentljsiiii^. In terms of 
the dimensionless filling n and the magnetization m de- 
fined in Eq. H2.6|l we may write 

fc| + fc^ = 2Tmla = ttN/L , (2.13a) 
k^-ki = 27rm/a = ttM/L . (2.13b) 

To study ferromagnetic symmetry breaking, we consider 
the change of the Gibbs potential due to a finite value 
of M . For convenience we introduce the dimensionless 
magnetization-dependent part of the Gibbs potential 



g{m) 



Nl 



r(n, to) — r(ri, TO = 0)] 



go (to) -I- /gi(m) + —g2{m) 



,(2.14) 



We suppress the dependence on the filling factor n, which 
is implicit in all quantities such as vp^ vq and /. In 
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exponential of Eq. (|2.1|l as follows, 

H - fiN - hM = Ko + V 

where 

ka 



and 



with 



i k(7 



£,ka ^ ek + Act - - crh . 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



In the language of many-body theory, A^- = ^a-ika, iO) is 
the (a priori unknown) self-energy due to the two-body 
interaction of our original hamiltonian 1)1. 1() for momenta 
k = ka- at the Fermi surface and for vanishing frequency. 
Following the usual procedurei^, the counterterms A^- 
can be determined order by order in perturbation theory 
by demanding that the self-energy generated by the sub- 
tracted interaction V in Eq. (|2.17l) vanishes for k — 
and uj = Q. Alternatively, the self-consistent determina- 
tion of the counterterms can also be implemented non- 
perturbatively within the framework of the renormaliza- 
tion groupi^ii^. 



B. Non-interacting limit 

Let us first consider the non- interacting limit [7 = 0, 
where Ag- = and 



no{^Ji,h) = -T^ln[l 



-/3(efc-Ai-(T?i)l 



(2.19) 



ka 



The expressions relating particle number and magneti- 
zation to chemical potential and magnetic field are then 



FIG. 2; (Color online) Density of states at the Fermi energy 
in units of for the dispersion 11.21 as a function of filling 
n for different values of 7 = it' /t, see Eas. llTT?71 ITTTI ITHl . 
The numbers close to the curves in a) give the corresponding 
values of 7. Note that for 7 > 1 the DOS remains finite for 
n — > ric from below, but diverges if n approaches ric from 
above. In contrast, for 7 < — 1 the DOS is finite if n ^ ric 
from above and is singular it n ^ ric from below. 

Eq. (|2.14|l we have normalized the Gibbs potential by 

the natural scale — = (—) of the kinetic energy. 
In the regime where the Fermi surface without symme- 
try breaking has only two pointsi^ it is then sufficient 
to set up the renormalized perturbation theory by intro- 
ducing two counterterms quadratic in the fermions with 
coefficients A| and A| and rewrite the operator in the 



N - ^ /(efc - /i - (T/i) , (2.20a) 

ka 

M ^ ^•^(^'^ ~^^-cTh) , (2.20b) 

ka 

where /(e) = [e^'^ -I- 1]^^ is the Fermi function. We de- 
note the corresponding Gibbs potential by ro(iV, M). In 
general it is not possible to calculate ro{N,M) analyti- 
cally. However, for weak ferromagnets, where |m| <C 1, 
we may expand F(A^, M) in powers of M. For the m- 
dependent part of the dimensionless Gibbs potential de- 
fined in Eq. (|2.14l) we obtain in the non-interacting limit 
at zero temperature, 

goim)=m^ + ^m'^ + 0{m'^) , (2.21) 
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FIG. 3: (Color online) Dimensionless coefficient C defined in 
Eq. 12.221 for the energy dispersion given in Eq. 11.211 in the 
regime — 1 < 7 < 1 and for arbitrary filling n. 



where the dimensionless coefficient C of the quartic term 
can be written in terms of the derivatives of the DOS at 
the Fermi energy^J^ 



1^0 



3 --^ 



(2.22) 



Here u'q and Vq are the first and the second derivative of 
the energy-dependent DOS v{£) defined in Eq. (|2.5|1 at 
^ = 0. Equation H2.22|l is valid in any dimension provided 
we use the D-dimensional DOS. Interestingly, if in ID the 
Fermi surface without magnetic field has only two points 
±fci?, then the constant C can be related to the cubic 
term in the expansion of the energy dispersion around 
the Fermi surface as follows 



C 



where C3 is defined by 



C3_ 



(2.23) 



vpq 



2W 



+ 'iq' + 0{q^) . (2.24) 
6 



The quadratic term in the expansion H2.24|l cancels on the 
right-hand side of Eq. (|2.22l) . For the dispersion H1.2|l we 
obtain in the regime where the Fermi surface consists of 
two points, 



, 1 + 47COs(7m) 



(2.25) 



1 -|- 7COs(7rn) 

A graph of this expression is shown in Fig. 13 For clarity, 
in Fig. 21 we also show some cuts for fixed 7 through the 
surface in Fig. O While for I7I < 1/4 the coefficient 
C is negative for all densities, for I7I > 1/4 there exists 
always a regime of densities where C is positive. In fact, 
by choosing n sufficiently close to ric we can fine tune C 
to assume any desired positive value. 
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FIG. 4: (Color online) Dimensionless coefficient C defined in 
Eq. 12.22^ for different values of 7 = 4t' /t as a function of the 
filling n. 



C. Hartree-Fock approximation 

To first order in U the self-energy generated by the 
subtracted interaction V is E(j — Un^^ — A^, which is 
independent of momentum k and frequency lu. Here 



Nl 



(2.26) 



From the requirement that the subtracted interac- 
tion does not generate a momentum- and frequency- 
independent self-energy, = 0, we find the countert- 
erms to first order in U, 



where we have defined 



(2.27) 



At -I- Ax Tjn]+ni 
A = — ' = U— = Un 



U 



2 



(2.28a) 
Um . (2.28b) 



The grand canonical potential is in this approximation 
given by 

fl{fi, h) « rioifi - A,h + A)- U{n^ - m^)NL , (2.29) 
and the corresponding Gibbs potential is 



(2.30) 



The function gi{m) defined in Eq. (|2.14() is therefore 
gi{m) = —rr? so that to first order in / we obtain 

g(m) = (l-/)m2 + :^m* + 0(m^/2) (2.31) 

In the regime where C > this leads to the usual Hartree- 
Fock scenarioiifi of a second order quantum phase transi- 
tion to a ferromagnetic state at /c = 1, as shown in Fig.O 
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FIG. 5: (Color online) Hartree-Fock approximation for the 
dimensionless Gibbs potential g{m) given in Eq. llOT^ for 
7 — 0.9 and n = 0.9, corresponding to C ~ 166. The numbers 
close to the curves are the corresponding values of the Stoner 
factor I. Note that for n — 0.9 the physical values for m are 
in the range |m| < 1 — n = 0.1. 




FIG. 6: Feynman diagram representing the correction to the 
Gibbs potential to second order in the on-site interaction. The 
wavy lines represent the bare interaction and solid arrows rep- 
resent the Hartree-Fock Green-functions with the indicated 
spin projections. 



Of course, extrapolating the Hartree-Fock approximation 
to J ~ 1 is an uncontrolled procedure, because in this in- 
teraction range there is no reason why the corrections of 
order P and higher in Eq. (|2.31|) should be small. In 
order to assess the validity of the Hartree-Fock scenario, 
we shall calculate in the following section the correction 
to the Gibbs potential to second order in /. 



III. GIBBS POTENTIAL TO SECOND ORDER 
IN THE INTERACTION 

The second order correction to the Gibbs potential is 
given by the Feynman diagram in Fig. |^ In real space 
and imaginary time this diagram represents the following 
expression. 



r2(iv, M) 



'Nl—- / dr 

^ J~l3/2 



where 



H^(x,r) = -G'ct(x,t)G'ct(-.t, -t) 



(3.1) 



(3.2) 



is the polarization bubble for spin a electrons without 
interactions. Here the real-space imaginary-time Green 
function is 



1 ^i{kx-LOiT) 



(3.3) 



where w; = 2tt{1 + i)T, I = 0, ±1, ±2, . . . are fermionic 
Matsubara frequencies. Carrying out the Matsubara sum 



and taking the zero temperature limit {(3 — > oo) and the 
infinite system limit [L — > oo) we obtain 



TT I a 



27r 



(3.4) 



with 

Gka{T) 



[e(r)e(a.) - e(-r)e(-a.)] 



(3.5) 

To make progress analytically, we linearize the energy 
dispersion within an interval —A < g < A around the 
Fermi points. Here A is an ultraviolet cutoff of the order 
of 1/a and we assume that 



m ^ — — A 

TT 



(3.6) 



which is necessary in order to justify the linearisation of 
the energy dispersion at the Fermi points. The dominant 
correction term to the Hartree-Fock approximation turns 
out to depend only logarithmically on A, so that it is not 
very sensitive to the numerical value of A. With these 
assumptions the integration in Eq. (|3.4|l can be carried 
out exactly. 



G.(.,r) « ^ 



1 



X + iv^T 



X — IVrrT 



(3.7) 
In this 



where Sr = sign(T) and v„ = d^k<j /dk\k=k^ 
approximation the polarization can be written as 

H,(a;,T) =nO(a:,r)+Hf'(a;,r)+H-2'=.(a;,r) , (3., 



6 



with the forward scattering contribution 



n°(x,r) = 



rRe 



(27r)2 

and the backscattering part 



(27r)2 a;2 + (uo-t)2 



(3.9) 



(3.10) 



Substituting Eq. (|3.7() into Eq. H3.1|l and taking the hmts 
f3 oo and L — > oo, we find 



2„3 



r2(A^,Af) = -NlU'cl 



(2^)4 



J-oo 



X Re 



l-2g|e*''^ + .gfe 



Re 



■ cos(2fc|x" 
' cos(2fc|X^ 



{x + TOi-r)2 
1 — 2g-f cos(Aa;) + g| 

1 — 251 cos(Aa;) + 

a;2 + {vitY 



{x + 



Re 



Re 



Ax I JlJliKx 



1 - 2.g^e'^^ + 
(a; + iv\TY 



+ cos{2k^x) cos{2kix) 



1 — 2gf cos(Aa;) + .g^ 
a;^ + (ftT)2 

I 



(a; + iv^r)^ 



1 - 2gj^ cos(Aa;) + g'f 
x^ + (wi'r)2 



(3.11) 



Here g„ — e^'""^'^ . The a::-integration can now be per- 
formed using the residue theorem. The resulting r- 
integration can then be carried out exactly. Keeping 
in mind that the difference v-^ — U| between the Fermi 
velocities is proportional to mU ^ we may approximate 
i)| vi K, vf 'va Eq. (|3.11|) . because the prefactor is 
already of order . In this approximation we obtain 

<?2(m) = ^ [r2(iV,A/) -r2(A^,0)] = X^f{m/\) , 



with 
/(^ 



(3.12) 

2|a;|)ln(l - 2|a;|) 



—x^ In |a;| — 3|x| — |a;|(l 
-(l-a;2)ln(l-|a:|) 

+ (1 + |a;| -2a;2)ln(l + 2|a;|) . (3.13) 

Because Eqs. (I3.12f) and (|3.13l) have been derived assum- 
ing that X ~ m/\is small (see Eq. (13.611 ) it is consistent 
to expand f{x) in powers of x, 



/(^) = 



-Inbl 



6|a:|3 -f -f 0(a;^) (3 -^4) 



Then we obtain the magnetization-dependent part of the 
normalized Gibbs potential to second order in / and to 
fourth order in m, 

511 2 3/2 ,3 

,. .4 

C 



gim) 



1 



13/2 



Inl A 

TO 



TO 



2A2 



^ + 0(|to|5,/3) 



(3.15) 



A reasonable choice for the ultraviolet cutoff A defining 
the interval where the linearization of the energy disper- 
sion is justified is A = n/a for n < 1/2 and A = {l — n)/a 
for n > 1/2. Our dimensionless cutoff A defined in 
Eq. H3.6|l is then A = vr^^minjn, 1 — n}. A graph of 
Eq. H3.15|l for this choice of A and n = 0.9 is shown in 
Fig. [71 The non- analytic terms proportional to In |to| 



and 



generated by the second order correction com- 



pletely change the Hartree-Fock scenario of a second or- 
der quantum phase transition depicted in Fig. |S1 From 
Fig. Ul it is obvious that the extrapolation of the second 
order correction to an interaction strength of the order 
of unity leads to a first order quantum phase transition, 
which for our choice of 7 = 0.9 and filling n = 0.9 oc- 
curs at a critical value /c w 0.6615. At this value of 
/ our function g{'m) develops three separate degenerate 
minima, which is the characteristic feature of a first or- 
der phase transition. The precise numerical value for the 
critical Ic depends on our particular choice of the cutoff 
A, so that the above value of Ic should not be taken too 
serious. However, with any reasonable choice of A the 
critical Ic is smaller than the Hartree-Fock result Ic = 1, 
suggesting that correlations can stabilize the ferromag- 
netic state in certain parameter regimes. For consistency, 
we should require that the magnetization to at the phase 
transition satisfies to < A ~ (1 — n)/7r, see Eq. (|3.6|l . 
From Fig. we see that this condition is only marginally 
satisfied, so that our approximations loose their quanti- 
tative accuracy once the curves in Fig. [7| develop minima 
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FIG. 7: (Color online) Gibbs potential g(m) to second order 
in the interaction given in Eq. 13.1511 with A = {l — n)/n. The 
values of the Stoner factor are written next to the curves. The 
parameters 7 = 0.9 and n — 0.9 as well as the scales are the 
same as in Fig. |3 



at finite m. 

It is instructive to examine the physical origin of 
the dominant logaritlimic correction proportional to 
In |to| in Eq. H3.15|l . This term arises from a product 
of two backscattering contributions to the polarization in 
Eq. which for small m contain also a component 

involving only small momentum transfers. In fact, if we 
are only interested in the non-analytic contributions to 
r2 {N, M) , we may replace Eq. H3.1|l by the simpler ex- 
pression 



-ising 



dT 



dx 



li]^' {x, r)n^'^'^ {x, t) + nf (x, T)nf ^ (x, t) 



2„3 



dT 



dx- 



is[2(A:| 



(3.16) 



where tq = {vf^)~^ and = {Ai:\m\v p / a)~^ . Using 
/c| — fcj = 27rm/a and expanding imder the integral sign 
for small m 

cos[47rmx/a] = 1 - -^^^m^ + 0{m^) , (3.17) 

it is easy to see from Eq. (|3.16f) by power counting that 
the coefficient of m? is decorated by a non-analytic cor- 
rection proportional to ln(Tm/To) = ln(aA/47r|m|). 



IV. DISCUSSION AND CONCLUSIONS 

In this work we have evaluated the magnetization- 
dependent part g{m) of the Gibbs potential V{N, M) 



to second order in the relevant dimensionless interac- 
tion strength / — Uvq. Our main result is an explicit 
expression for g(r7i) up to order given in Eq. H3.15|l . 
which contains non-analytic corrections proportional to 



'In I 



and 



When extrapolated to interactions 



/ of the order of unity, these corrections imply that the 
paramagnetic-ferromagnetic quantum phase transition in 
ID, if it exists, must be first order. Of course, the ex- 
trapolation of the weak coupling expansion to values of 
/ of the order of unity is an uncontrolled procedure, so 
that with our method we cannot proof the existence or 
the absence of a ferromagnetic ground state in ID. How- 
ever, the numerical density-matrix renormalization group 
calculations by Daul and Noacls^ suggest that in a cer- 
tain regime of 7 = At' /t and fillings n a ferromagnetic 
ground state indeed exists. On the other hand, these au- 
thors found numerical evidence that the paramagnetic- 
ferromagnetic quantum phase transition in ID is second 
order, which is not supported by our calculation. 

Given the fact that the metallic state in ID is a Lut- 
tinger liquid, one could have expected that the pertur- 
bative expansion of the Gibbs potential contains non- 
analytic corrections. Surprisingly, according to a recent 
calculation by Betouras, Efremov, and Chubukov^^ sim- 
ilar non-analytic corrections exist even in higher dimen- 
sions. These authors evaluated the magnetic-field de- 
pendent part of the interaction-correction to the grand 
canonical potential to second order in U in three and 
two dimensions. They found that the susceptibility is 
proportional to ft,'^ln(|/i|) in three dimensions, and pro- 
portional to \h\ in two dimensions. These corrections im- 
ply non-analyticities in the corresponding magnetization- 
dependent part of the Gibbs potential, leading to the 
breakdown of the Hartree-Fock approximation in dimen- 
sions D > 1, in agreement with the work of Belitz and 
Kirkpatrick^. 

Our result (|3.15|) for the Gibbs potential implies that 
the magnetic susceptibility at weak coupling vanishes for 

TO — > as 



x(to) 



d'^g{m) 



dr 



2vo 



P ln(A/|TO|) 



(4.1) 



Hence, in the paramagnetic state the zero-field suscep- 
tibility x(™ = 0) vanishes for any finite value of the 
interaction, while it approaches the finite value vq in the 
non-interacting limit. In one-dimensional metals other 
quantities are known to exhibit a similar discontinuity. 
For example, the density of states of a Luttinger liquid 
vanishes at the Fermi energy, while in the absence of in- 
teractions it is finite-'^^. We conjecture that higher order 
corrections neglected in Eq. (|3.15() will transform the log- 
arithmic singularity in Eq. (|4.1|l into a power law with 
interaction-dependent exponent. 

It is important to emphasize that a ferromagnetic 
ground state in a ID lattice model with hopping be- 
yond the nearest neighbors does not contradict the Lieb- 
Mattis theorem-"^. Ferromagnetic ground states of the 
t — t' Hubbard model can therefore not be ruled out a 
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priori. Whether such a ferromagnetic ground state is 
relevant to explain recent measurements of the conduc- 
tance anomaly in quantum wiresiS^ is not clear at this 
point. Very recently Klironomos et alP'^ pointed out that 
electron-electron interactions in quantum wires induce 
deviations from the strictly one-dimensional geometry, 
in which case the Lieb-Mattis theorem does not apply 
and a ferromagnetic ground state is in principle possible. 

A detailed evaluation of the second order correction 
to the Gibbs potential in the regime where the Fermi 
surface without interaction consists of four points still 
remains to be done. This calculation is more compli- 
cated than in the case of two Fermi points, because a 
self-consistent treatment within renormalized perturba- 
tion theory requires the introduction of four countert- 
erms, one for each Fermi poinfei^*i^. VoUhardt et at— 
pointed out that a large asymmetry in the DOS with a 
peak at the lower band edge tends to favor ferromag- 
netism by minimizing the increase of kinetic energy due 
to the spin polarization. Given the fact that in our model 



for I7I > 1 the DOS exhibits a one-sided singularity at 
the critical filling ndj) where the number of Fermi points 
changes discontinuously (see Eqs. (|2.81 12.9(1 and Fig. [SJ, 
we suspect that for fillings n close to ric a ferromagnetic 
ground state can be stabilized even for values of the in- 
teraction strength / that are substantially smaller than 
unity. Whether this hypothesis is correct requires a de- 
tailed calculation, which is beyond the scope of this work. 
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